Polymer-induced tubulation in lipid vesicles 
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A mechanism of extraction of tubular membranes from a lipid vesicle is presented. A concentration 
gradient of anchoring amphiphilic polymers generates tubes from bud-like vesicle protrusions. We 
explain this mechanism in the framework of the Canham-Helfrich model. The energy profile is 
analytically calculated and a tube with a fixed length, corresponding to an energy minimum, is 
obtained in a certain regime of parameters. Further, using a phase-field model, we corroborate 
these results numerically. We obtain the growth of tubes when a polymer source is added, and the 
bud-like shape after removal of the polymer source, in accordance with recent experimental results. 

PACS numbers: 87.16.ad,87.16.Wd 



Introduction.- Vesicle shape transformations, like 
tubulation, play an essential role in cellular transport. 
Both energy and matter are transferred continuously 
throughout the membrane, a lipid bilayer with proteins 
and other macromolecules anchored on it, which sur- 
rounds the cell and most of its internal organelles [l|. 
As part of cellular dynamic processes, these membranes 
adopt different shapes in order to exchange matter with 
their surroundings. Many possibilities appear here, from 
budding and eventual fission of small transport vesicles 
to formation of large tethers connecting distant or- 
ganelles, as in the Golgi apparatus and the endoplasmic 
reticulum [3], or even between different cells ^]. 

The formation of these tethers can be driven by ap- 
plying a point-like force to the membrane [a, |6[. Un- 
derstanding the nature of this force is of major impor- 
tance. There are different mechanisms leading to such 
a tubulation phenomenon, as for instance the growing 
of microtubules pushing the membrane from inside _7|, 
and the extrusion due to a hydrodynamic flow 8] . Other 
works have experimentally studied the force generated by 
molecular motors pulling membrane tubes in vitro [9| and 
by optical tweezers [l3| ■ Tsafrir et al. studied the tubula- 
tion induced in highly oblate vesicles by the anchoring of 
amphiphilic polymers |ll| without any directed force. In 
those experiments, macromolecules containing hydropho- 
bic groups were administered in the surroundings of a gi- 
ant oblate vesicle. Those molecules diffused in the bulk 
and eventually anchored the membrane inducing a local 
spontaneous curvature by the mechanism of hydrophobic 
insertion [l2| , leading to the formation of one or several 
buds. Even, under certain circumstances, those buds can 
grow into long tubular structures (see Fig. [T]). 

Here, motivated by those experimental results [llj , we 
present a theoretical treatment of a novel mechanism of 
tube extraction, due to the generation of spontaneous 
curvature by anchoring macromolecules which are dis- 



tributed along a gradient of concentration, maintained 
by a source. We show that the elongation of a bud into 
a tube may be, under certain circumstances, energeti- 
cally favorable when a polymer concentration gradient is 
present. We obtain that these tubes do not grow indef- 
initely, but stop at a certain length. In the first part of 
this Letter, we analyze the problem in the framework of 
the Canham-Helfrich model [IJ, [ij] with a simple geome- 
try and a stationary linear concentration profile, in order 
to solve analytically and to understand qualitatively the 
tube formation. In the secondpart, we use a phase-field 
model for the bending energy 15| , coupled with a station- 
ary polymer concentration profile, to study the problem 
numerically. 
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FIG. 1: Tube formation in a highly oblate vesicle after local 
addition of multianchor polymer. Experimental images from 
Ref. [11. 



Theoretical treatment.- Let us assume that a hemi- 
spherical bud of radius R is already formed out from a gi- 
ant vesicle. This bud has not to be necessarily stable, ap- 
pearing for instance just as a vesicle rim fluctuation [llj. 
The mother vesicle is large enough to be considered as a 
lipid reservoir during the whole extraction process, so we 
can study the isolated tube on its own. In addition, there 
are amphiphilic macromolecules (e.g. polymers as in (llj) 
in the volume outside the vesicle. We assume here, in or- 
der to find analytical estimations for the extraction of 
membrane tubes due to an inhomogencous polymer con- 
centration in the bulk, that these macromolecules are ap- 
plied in a line located at a distance Zp from the mother 
vesicle (see Fig. [2]) and follow a linear stationary pro- 
file for this concentration, p(r,t) = po {1 — \z/zp — 1|), 
where Tappi = (0,0, Zp). 

The bending energy of a membrane treated as an elas- 
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FIG. 2: Geometrical sketch of the system: a cyhndrical tube 
of radius R and length L, with a hemispherical cap of radius 
R. A linear polymer concentration gradient is p(z) is also 
outlined in the figure (dashed line). 



tic sheet is given by the Canham-Helfrich Hamiltonian 
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where ci and C2 are the principal curvatures of the sur- 
face, cq the spontaneous curvature and k the bending 
modulus. Since the tube is connected to a lipid reser- 
voir, we have also included a surface tension term, where 
(T is the surface tension of the membrane. A pressure-hke 
term could also be introduced, but its effects are sublead- 
ing [5|. Using a cylindrical tube with a hemispherical cap, 
as sketched in Fig. [21 Eq. [T] can be written as 
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where Sl and Sr denote the cylindrical and the henri- 
spherical parts of the tube (Fig. [2]) with areas Al and A^, 
respectively. The spontaneous curvature of the mem- 
brane is coupled with the local concentration since the 
hydrophobic anchor groups of the polymer tend to insert 
themselves into the bilayer, acting thus as a wedge \W{ . 
A linear coupling, as stated previously in a similar sys- 
tem [1^ [l3|, has shown to be justified, the spontaneous 
curvature field outside the vesicle being co(z) = cq p{z). 
Therefore, in the case where the polymer gradient in 
space is set to be linear, we get for the tube energy Eq. ^ 
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where R — 1 sets the length scale, a — a/ k,, and we de- 
fined £, — Cq pq/ Zp, as the slope of the linear spontaneous 
curvature profile. 

In Fig. [3] (inset) we show how the tube energy Eq. ^ 
looks like as a function of the length L of the tube for 



different slopes of the spontaneous curvature profile, ^, 
when the tension a is negligible. Note that the energy 
Eq. ([3]) is cubic in the tube length. The energy extremes 
correspond to two equilibrium lengths: one stable length 
corresponding to an energy minimum, L* = 1/^ —1-1- 
\/'2/£, — 7r/2 + 1; and another smaller length being unsta- 
ble, Lc = 1/C- 1 - v/2/C- 71-/2 + 1. The larger the slope 
of the linear concentration profile is, the smaller these 
lengths are. For ( > ^"'^ = 4/7r(l - ^/l - 7r/8) - 0.28 
(see Fig. |3|) , theres is only a local minimum and no local 
maximum of the energy for positive lengths (Fig. [3]). This 
means that the initial bud becomes unstable against the 
formation of a tube. For smoother concentration gradi- 
ents (^ < C'^"*), the initial bud, in order to grow up to 
its preferred length L*, needs to cross an energy barrier 
AEq, or to start with a certain initial length larger than 
Lc (see Fig.|3|). For very steep slopes of the concentration 
gradient f > ^"'^'^ = 4/7r(l + y/l-n/8) ~ 2.27, the lo- 
cal energy minimum at L* disappear, and no stable tubes 
can be found (see Fig. 2]). 
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FIG. 3: Energy vs. tube length for different values of the 
rate of polymer molecules added to the system, i.e. the slope 
of the linear polymer profile ^ (inset). Initial buds need to 
overcome an energy barrier AEq to elongate into tubes for a 
certain range of values of ^. Then, the elongated tube needs 
to overcome another barrier e in order to be reabsorved by 
the mother vesicle. The surface tension here is negligible. 

For an initially formed bud-like fluctuation {Lmi — 
in Fig. [4]), when we increase the slope of the sponta- 
neous curvature, ^, the bud cannot grow unless the criti- 
cal slope, £,^"*^, is reached, and then a finite- length tube is 
formed out (dashed line in Fig. 2]) . This transition is dis- 
continuous in the value of the stable length of the tube. 
In addition, we see that the tube length is maximum for 
f'^"*, decreasing for steeper profiles, <^ > ^"^"^ Eventu- 
ally, at ^'"^'^ the equilibrium tube length vanishes, and 
buds become again unstable. This transition is therefore 
continuous. 

For non-vanishing but relatively small values of the 
surface tension, a, the stable tube length decreases as 
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FIG. 4: Stability diagram of tube formation as a function of 
the initial tube of length Lini (iini = corresponds to a bud), 
and the dimensionless slope of the spontaneous curvature gra- 
dient 5. The dashed line correspond to the stable tube length, 
L* , for a given concentration gradient. 



L*{a) = i*(0) -ct/(^2 ^i _ ^/2 + 2/0), where L*{0) is 
the tube length for zero tension. For tensions larger than 
a critical tension CTc = 4 — 7r/2, bud to tube transition 
dissapear, altough stable tubes may be formed out from 
a tube with a finite given length. 

Numerical model.- In the experiments by Tsafrir et al. 
[llj , polymer molecules diffuse from a source and eventu- 
ally anchor the membrane, inducing a local spontaneous 
curvature which modifies the vesicle shape. In order to 
study dynamically this process, we present here a numer- 
ical model dealing with a stationary polymer profile and 
the Canham-Helfrich energy without any assumption on 
the tube geometry. A systematic study of other station- 
ary and non-stationary concentration profiles, and the 
diffusion of the polymer molecules and their anchorage 
on the membrane, is beyond the scope of this Letter, 
and will be presented elsewhere llSil . 

We use a phase-field method [3, HI, M, EH, M, El 
to numerically deal with the bending energy of the 
membrane. Within this approach, the Canham-Helfrich 
Hamiltonian Eq. ^ is written as a dynamic function of 
a field, (/i, whose level-set {x : 0(x) = 0} locates the 
membrane position at each time. The dynamic equation 
for the phase-field, and hence for the membrane shape, 

is [m 



-£ = 2V^[{3q)^-l + 2eCoip{x))cj,)^(b,p] 



- e^V^<P[(l),p]+e^aV^(j)y 
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where ^[(j){x),p{x)] 



l)i<l>-eCoip{x))) 



(x), Cq{p{x)) — cop(x) is the spontaneous curva- 
ture induced by the local concentration p(x) of polymer 
molecules anchored on the membrane, and e is a small 
parameter related to the width of the interface. This dy- 
namic equation includes a tube surface tension a, and 
conserves locally the inner volume, due to the use of a 



relaxational model-B-like conserved dynamics. 

The assumption we did before of a linear stationary 
polymer profile is relaxed at this point. Since our aim 
in this Letter is to show how an inhomogeneous sta- 
tionary polymer concentration profile is a possible mech- 
anism of tube formation, we are going to use in our 
numerical treatment a Gaussian stationary profile such 

as p{x) — po/S\/27rexp ( — jx — Xappil /2i;^j, where 
S is the standard deviation, related with the width of 
the distribution, and po is the total number of poly- 
mer molecules introduced in the system. Our problem is 
then reduced to numerically solve Eq. ^ using this sta- 
tionary polymer concentration for different sets of initial 
conditions and parameters. We used a standard finite- 
difference scheme for the spatial discretization and an 
Euler method for the time-derivatives |13| . In all the re- 
sults shown in this Letter we used the value of the small 
parameter, e, equal to the mesh size of the lattice. The re- 
sults are robust under variations of this parameter. Both 
the time and space discretizations are chosen in order 
to satisfy the Courant-Friedrichs-Lewy stability criterion 

The parameters which are relevant in these simulations 
in order to study the growth and the shrinkage of tubes 
are: the length of the initial tube, and the characteristics 
of the polymer gradient profile. We fixed the position 
of the polymer source to be Zp = 7.5 i?, where R is the 
radius of the initial tube formed by a cylinder and a hemi- 
spherical cap. Due to the axisymmetry of the problem, 
the integration is performed in a two-dimensional lattice, 
whose size is 80 x 20 throughout this Letter. 

Num,erical results and discussion.- In Fig. [T] we show 
three snapshots of the experimental results from Ref. [ll| , 
where a fluctuating oblate vesicle undergoes a shape in- 
stability which forms an initial bud and eventually grows 
into a tube, due to the anchorage of amphiphilic poly- 
mer molecules. The presence of the polymer molecules 
enhance fluctuation of the vesicle rim (Fig. [Ha)). One 
can divide the dynamics of this process in four differ- 
ent regimes. The first one is the formation of the initial 
bud and suppression of rim fluctuations (Fig. [Ijb)). The 
second regime is when the tube starts to grow up to a 
certain length (Fig.[IJc)). Then, in the third regime the 
polymer source is shut down (Fig. [5fa)) and tubes shrink 
reaching a metastable bud- like shape (Fig. [SKb)). In the 
last regime, they eventually disappear reabsorbed by the 
mother vesicle. Within our model, we find, as shown in 
Fig-EJc), the shape of a large tube, which is qualitatively 
in good agreement with the experimental results. For a 
large range of values of the standard deviation, S, the 
tube length we get is essentially the same. In other words, 
when the source of polymer molecules is shut down, the 
concentration profile gets stretched as time goes by, but 
the length of the stable tube continues being the same. 

Then, after some time, the polymer molecules are more 
homogeneously distributed, and a new configuration of 
short length is found (Fig. [5Kd)). These buds are also in 
agreement with those found by Tsfarir et al. |II| , and we 



can quantitatively compare them by measuring the ratio 
between their width and length to be around 0.3. 
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FIG. 5: Tubes extruded from a vesicle by a non-homogeneous 
polymer profile. Comparison between the experimental re- 
sults from Ref. [lj| (a,b) and the phase-field integrations (c,d). 
For short times (small standard deviation E) long tubes are 
obtained (a,c), and for long time (wider distributions) buds 
appear (b,d). The resulting profile for the spontaneous cur- 
vature is shown for the phase-field integrations. We choose 
C = 0.3, S = S.5R (c), and C = 0.3, E = 107? (d). 

In other recent experiments, Roux et al. [25J gener- 
ated dynamin-coated membrane tubes. They observed 
the growth of these tubes by the addition of this protein. 
Further addition of GTP molecules to these tubes induces 
fission by a conformational change of dynamins. Altough 
our model is capable of explaining the initial tubulation 



by the generation of membrane curvature by dynamin, it 
does not consider topological changes such as fission. An 
extension of our model including the Gaussian curvature 
term in the Canham-Helfrich Hamiltonian eq. ([T]) will be 
published elsewhere [2(a] . 

Conclusions.- In this Letter we showed how a nonho- 
mogeneous concentration of amphiphilic molecules leads 
to tube formation from vesicles. The mechanism is ex- 
plained in the framework of the Canham-Helfrich model 
for the membrane energy, where the macromolecules an- 
chored on the membrane induce a local spontaneous cur- 
vature. First, we analyzed the energetics of such a sit- 
uation with a simplified geometrical scheme and a lin- 
ear profile for the concentration, and showed analytically 
under which conditions a tube may be formed from a 
vesicle. Afterwards, we used a phase-field model to nu- 
merically calculate the shape of tubes under no geometric 
assumptions with a Gaussian profile for the polymer con- 
centration. These numerical treatment led us to find the 
tube shape, in agreement with those found experimen- 
tally by Tsafrir et al. [11|. Two qualitatively different 
stable shapes are found depending on the characteris- 
tics of the polymer distribution: long tubes for initially 
narrow polymer distributions, and short buds for wider 
ones. 
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